Graviton mediated photon-photon scattering in general relativity 
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In this paper we consider photon-photon scattering due to self-induced gravitational perturbations 
on a Minkowski background. We focus on four- wave interaction between plane waves with weakly 
space and time dependent amplitudes, since interaction involving a fewer number of waves is excluded 
by energy-momentum conservation. The Einstein-Maxwell system is solved perturbatively to third 
order in the field amplitudes and the coupling coefficients are found for arbitrary polarizations in the 
center of mass system. Comparisons with calculations based on quantum field theoretical methods 
are made, and the small discrepances are explained. 
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I. INTRODUCTION 



As is wellknown, photon-photon scattering can occur 
due to the exchange of virtual electron-positron pairs, 
as described by QED or modifications thereof, see e.g. 
yj, and may even lead to collective photon phenomena 
2] . Photon interactions via the quantum vacuum, some- 
times involving deviations from the standard model, has 
recently been much discussed in the literature due to 
advances in experimental technologies (see e.g. |3|) as 
well as new theoretical insights ,4]. Moreover, photons 
also interact gravitationally, although this effect has been 
much less studied. Purely general relativistic treatments 
of electromagnetic wave interactions have been made re- 
sulting in exact solutions, see e.g. but these calcu- 
lations are very different from the pure scattering pro- 
cesses, and do not address the interaction at the single 
photon level. On the other hand, it is not clear to what 
extent calculations of the gravitational cross-section us- 
ing quantum field theoretical methods 00] (see also @) 
are consistent with classical general relativity. In order 
to shed light on this issue, we will consider the interac- 
tion of four electromagnetic (EM) waves on a Minkowski 
background, which is the lowest order scattering pro- 
cess consistent with energy-momentum conservation. By 
studying the classical Einstein-Maxwell system, but ig- 
noring terms that do not correspond to pure scattering 
(e.g. frequency shift terms) we will attempt to make con- 
tact between the classical and quantum field theoretical 
picture. Calculating the classical coupling coefficients be- 
tween waves of different polarizations, corresponding to 
the scattering amplitudes in quantum field theory, we are 
able to compare the classical cross-section with that of 
quantum field theory 0,0- While the results are approx- 
imately equal for small scattering angles 8, we find that 
there are significant differences for large 9. The likely 
source behind this discrepancy is that the quantum field 
theoretical calculation Q used the matrix scattering am- 



plitude in order to define the interaction potential. As 
shown by Ref. 0, such a procedure is not able to fully 
reproduce the general relativistic potential. 

Finally we note that while gravitational photon-photon 
scattering is weaker than the QED scattering in most 
cases of physical interest 0, Q , it should be noted that 
in the long wavelength limit, actually the gravitational 
cross-section is larger than that due to QED. 



II. THEORY AND RESULTS 

We employ units such that the speed of light and 
Planck's constant are c = H = 1, and use metric sig- 
nature (— , +, +, +). Tetrad indices a, 6, . . . run from to 
3 and a, /3, . . . from 1 to 3. Coordinate indices fi, v, ... go 
from to 3. 

Assuming plane waves and denoting the interacting 
waves A, B, C and D the total electric field is given 
by 



E 



(a^ett-M** + c.c. 



A,B,C,D, 



*Also at: Centre for Fundamental Physics, Rutherford Appleton 
Laboratory, Chilton, Didcot, Oxon OX11 OQX, U.K. 



where c.c. denotes the complex conjugate. A similar ex- 
pression holds for the magnetic field. Moreover we as- 
sume that the amplitudes have a weak dependence on 
space and time, i.e. \d^E(x u )\ <C \k^\ \E(x u )\. The pres- 
ence of the EM fields induce a perturbation, h^ u , of the 
flat background metric, f^„, enabling energy exchange 
between the modes. A generic frame, orthonormal to 
quadratic order in the field amplitudes (linear order in 
the metric perturbation), is chosen as 



-/ioi, 1 — -^hn, — -hi2, — -his 

2^02, —2^12, 1 — 7:^22, — 2~fr-23 



(1) 
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where e a = e^9 M . The matching condition corresponding 
to energy and momentum conservation is given by 



kg 



k c 



(2) 



In the center of mass system all frequencies are equal and 
the waves are counterpropagating pairwise, i.e. the wave 
vectors satisfy \la = — kg, kc = — ku. We make the 
following ansatz for the metric perturbations 



2lL0t 



h a p = h a p(x,z,t)e lu + h* af3 (x, z, t) e 
+h a/3 (t, x, z) e-^A+kcy* + ~ h * a/} (t) X) z) e4 (k^+k c) .x 
+h al3 (t, x, z) e -^A-k c y* + ^ fa x> z ) e >(kA-k c )-^ 

This is not the most general ansatz, however it is suffi- 
cient to give all terms corresponding to resonant energy 
exchange between the modes. All the components h a @, 
h* a p-> h a p, h* a p, h a p, h* a g, except twelve of them, are de- 
termined by the field equations, G a b = KT„t, where G a b is 
the Einstein tensor and T a b the energy-momentum ten- 
sor. In the present case the only contribution to the 
energy-momentum tensor is given by 

Tab — F^Fbc — \g a bF cd F c d . 

The undetermined coefficients in the metric ansatz are 
set to zero using the generalized Lorentz condition. Note 
that the non-zero coefficients are of quadratic order in the 
field amplitudes. Using Maxwell's equations V[ a F bc j = 
and V a F ab — j b , where F ab is the electromagnetic field 
tensor and j° the four-current density, we can derive the 
following wave equations 



DE a 



t O (30 e a i3 7 

e a ^C^ e a E 



J 0i 



(3) 



(4) 



Here the wave operator □ = eo-eo+V- V, which coincides 
with the D'Alembertian operator in Euclidian space, and 
C c ah are commutation functions for the frame vectors sat- 



isfying 



Cab e c- 



5e, Jb, Pe and ps are the ef- 



fective currents and charges due to the inclusion of the 
gravitational field given by 



3B 

Pe 
Pb 



e { L 00 1 



01 



B s ) 



^O0B° 



1 0a D 



a/3 7 -p0 p 
e L afi^l 



(5) 

(6) 
(7) 
(8) 



where T% c are the Ricci rotation coefficients. The effec- 
tive currents and charges will be of cubic order in the 



field amplitudes. Eliminating the magnetic field from 
the wave equation @ by using Faraday's law to lead- 
ing order, k x E = wB, and neglecting terms of order 
four or higher in the field amplitudes will result in three 
categories of terms. 

1 . Nonresonant terms that will vanish after averaging 
over several wavelengths and time periods. Inter- 
action due to these terms are not consistent with 
energy momentum conservation. Most terms be- 
long in this category. 

2. Phase shift terms which are resonant but give rise 
to phase shifts rather than scattering. These terms 
typically contain a certain amplitude together with 
its complex conjugate, i.e. terms of the form 
EaE* a Eb- 

3. Resonant scattering terms containing all three wave 
amplitudes according to the energy-momentum 
conservation condition (J2J. 

Based on the classification of terms above we restrict 
ourselves to include only the resonant scattering terms 
and introduce polarization states perpendicular to the 
wave vectors as shown in figure 1. Note that the E 2 and 
the E x directions coincide. We thus have 

E A1 = - cos {9/2) E A+ , E A3 = sin (6/2) E A+ , 

E B i = cos (9/2) E B+ , E B3 = - sin (9/2) E B+ , 

Eci = - cos (9/2) E c +, E C3 = - sin (9/2) E c +, 

E m = cos (9/2) E D+ , E D3 = sin (9/2) E D+1 

with E A2 = E Ax , E B2 = E Bx , Eci = Ecx , E B2 = E Bx , 
and 

k Ax = sin (9/2) u>, k Az = cos (9/2) ui, 
k Bx = — sin (9/2) ui, k Bz = — cos (9/2) to, 
kcx = -sin (0/2) w, k Cz = cos (9/2) u>, 
k Bx — sin (0/2) u, k Bz — — cos (0/2) w. 

After some lengthy but straightforward algebra we end 
up with the following coupling equations describing the 
evolution of the wave amplitudes to leading order 

OE A+ = FiE* B+ Ec+E B + + F 2 E B+ EcxE B x 

+F 3 E Bx EcxE B + + F i E* Bx Ec+E Bxl (9) 

DE B+ — FiE* A+ Ec+E B+ + F 2 E* A+ Ec x E Bx 

+F 3 E* Ax Ec+E Bx + F i E* Ax Ec x E B +^0) 

DEc+ — FiE* D+ E A+ E B+ + F 2 E* D+ E Ax E Bx 

+F 3 E D . x E Ax E B+ + F±E* Dyc E A+ E Bx (11) 

DE B+ = FiE c+ E A+ E B+ + F 2 E c+ E Ax E Bx 

+F 3 E Cx E A+ E Bx + F4E Cx E Ax E B+ (12) 
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FIG. 1: Polarization directions and the definition of the scat- 
tering angle 9. 



FIG. 2: The classical(solid line)and quantum field theoretical 
(dashed line) cross-sections in arbitrary units as a function of 
scattering angle. 



where □ = d /dt — d /dx — d /dz and 



Fi = 

F 3 = 
Fi = 



K ( 3 + cos 2 



1 — COS 



4k (2 + cos 2 9 + cos ( 

1 + cos 9 
4k (2 + cos 2 9 - cos ( 

1 — cos 9 



(13) 



For symmetry reasons OEax , DEbx , OEcx and OEdx 
can be found from (j9" ]) - H12|) respectively by interchanging 
+ and x . The coupling coefficients only depend on the 
scattering angle 9, and in the limit when 9^0 both 
Fi and F 4 become infinite while F 2 — > — F 3 . The small 
angle divergence in Fx and F 4 is a consquence of the 
infinite range of the gravitational force. However, the 
coefficients F 2 and F3 must remain finite for all angles, 
as those coefficients not only describe scattering an angle 
9, but also correspond to a change in the polarization 
state. 

In order to check the consistency of our results, we 
assume long pulses, i.e. □ w —2iu>dt, such that the time 
derivative of the total energy density, e to t = J2 n (\Fn+\ 2 + 
\E nx | 2 ) (where the sum is over A, B, C, D), can be easily 
calculated. Carrying out the sum, it is found that all 
the scattering terms cancel, and thus we deduce that the 
evolution equations are energy conserving. 

Next we rewrite (|9"j)- (|12[l in terms of the vector poten- 
tials, which rescales the coupling coefficients (|13fl by a 
factor u) 2 . Noting that the rescaled coupling coefficients 
corresponds to the scattering amplitudes, and following 
Ref. [13, we find that the unpolarized differential cross- 
section can be calculated as 



da_ 

an 



128w 2 (2wY 



(14) 



where the square of the scattering matrix amplitude av- 
eraged over all polarization states is given by 



sin 4 9 



8 9 + 28 cos 6 61 + 70 cos 4 9 

-28 cos 2 9 + 129) . (15) 



This result should be compared with the same quantity 
calculated by quantum field theoretical methods, i.e. Eq. 
(15) in Ref. Q- It turns out that the differential cross- 
sections agree in the limit 9 <C vr/2, but as seen from 
Fig 2, where the classical and quantum field theoretical 
expressions are shown, the two expressions differ slightly 
in general. In order to resolve the difference more accu- 
rately, the ratio of the cross-sections are shown in Fig. 
3, where one should note the agreement for small an- 
gles. However, for general angles the expressions clearly 
disagree, and it is natural to ask what causes this dis- 
crepancy. To answer this question we note that Ref. 
has used the matrix scattering amplitude to determine 
the interaction potential. As demonstrated by Ref. 0, 
however, such a procedure is not sufficient to fully re- 
produce the general relativistic potential. As the general 
relativistic deviation from Newtonian behavior becomes 
more pronounced for large scattering angles, this explains 
the deviation for general angles, but also the agreement 
in the small angle limit. 



III. SUMMARY AND DISCUSSION 

Comparing the calculated cross-section for gravita- 
tional photon-photon scattering Eq. I|14|) with that from 
QED photon-photon scattering (due to exchange of vir- 
tual electron-positron pairs), we find that they have dif- 
ferent frequency dependence. The former is proportional 
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FIG. 3: The ratio of the classical and quantum field theo- 
retical cross-sections <n and 02, respectively, as a function of 
scattering angle. 

to u! 2 while the latter is proportional to oj 6 10]. Not- 
ing from Eq. (JUJ that da/dn ps 16(2n) 2 Lp\ 2 (letting 

sin 9 ~ cos ~ 1/ \/2) , where L p is the Planck length and 
A the wavelength, and comparing with the QED expres- 
sion for the cross-section (e.g. Ref. IToj). we find that 
the QED and gravitational cross-sections are comparable 
for frequencies 

c~10 3 C (L>A?) 1/2 , (16) 

where ro is the classical electron radius, A c is the Comp- 
ton wavelength, and we have reinstated the speed of 
light c. Thus the gravitational effects become the dom- 
inant contribution to the cross-section for frequencies 
uj ~ 30rad/s and lower. Still, the cross-section is very 



small, and we need extremely large photon densities for 
gravitational photon-photon scattering to influence the 
dynamics. Situations that could be of interest to study 
in more detail involve the dense photon gas surrounding 
pulsa rs 1 11! . as well as the photon gas in the early uni- 
verse • Furthermore, we note that if the energy densi- 
ties are sufficiently high 13], the timescales for nonlinear 
evolution will not be determined by the cross-section, 
even if the spectrum is strongly incoherent. Instead the 
characteristic time-scale must be found from weak tur- 
bulence theories Using the so called random phase 
approximation, the phase dependence can be integrated 
out, and evolutio n eq uations for the spectral energy den- 
sities are derived [l4l ]. 

Graviton mediated photon-photon scattering share 
many parameter similarities with photon-graviton pair 
conversion [15]. While it is possible that gravitational 
photon-photon scattering may have applications to astro- 
physics and/or cosmology, the effect can typically be ne- 
glected compared to other effects, such as QED photon- 
photon scattering or, in the presence of matter, interac- 
tion with charged particles. Thus our main aim here has 
been to make an explicit comparison with the quantum 
field theoretical result (see Fig. 2), which show a slight 
deviation from our general relativistic cross-section. As 
seen in Fig. 3, the deviation vanishes in the limit of small 
scattering angles. We trace the difference between the 
quantum field theoretical and the classical result to the 
difficulty in determining the general relativistic interac- 
tion potential from the matrix scattering amplitude, as 
done by Ref. pfl. An interesting problem, which is a 
project for further research, is to investigate whether a 
quantum field theoretical calculation can be improved to 
incorporate a fully general relativistic interaction poten- 
tial. 



[1] W. Heisenberg H. and Euler H., Z. Physik 98, 714 (1936); 
J. Schwinger J., Phys. Rev. 82, 664 (1951). G. Brodin, 
M. Marklund and L. Stenflo, Phys. Rev. Lett. 87, 171801 
(2001). 

[2] M. Marklund, G. Brodin, and L. Stenflo, Phys. Rev. Lett. 
91, 163601 (2003); M. Marklund, G. Brodin, L. Stenflo, 
and P. K. Shukla, Phys.Scripta T107, 239 (2004); M. 
Marklund and P. K. Shukla, Rev. Mod. Phys. 78, 591 
(2006). 

[3] G. A. Mourou, T. Tajima, and S. V. Bulanov, Rev. Mod. 
Phys. 78, 309 (2006); Y. I. Salamin, S. X. Hu, K. Z. 
Hatsagortsyan, and C. H. Keitel, Phys. Rep. 427, 41 
(2006). 

[4] E. Lundstrom et al., Phys. Rev. Lett. 96, 083602 (2006); 
E. Zavattini et al., Phys. Rev. Lett. 96, 110406 (2006); 
R. Rabadan, A. Ringwald, and K. Sigurdson, Phys. 
Rev. Lett. 96, 110407 (2006); B. A. van Tiggelen, G. 
L. J. A. Rikken, and V. Krstic, Phys. Rev. Lett. 96, 
130402 (2006); D. B. Blaschke, A. V. Prozorkevich, C. 



D. Roberts, S. M. Schmidt, and S. A. Smolyansky, Phys. 
Rev. Lett. 96, 140402 (2006); A. Di Piazza, K. Z. Hat- 
sagortsyan, and C. H. Keitel, Phys. Rev. Lett. 97, 083603 
(2006); J. T. Mendonga, J. Dias de Deus, and P. C. Fer- 
reira, Phys. Rev. Lett. 97, 100403 (2006); R. Schiitzhold, 
G. Schaller and D. Habs, Phys. Rev. Lett. 97, 121302 
(2006). 

[5] J. B Griffiths, Colliding Plane Waves in General Rela- 
tivity, Clarendon Press, Oxford (1991) 

[6] B. M. Barker, S. N. Gupta and R. D. Haracz, Phys. Rev., 
149, 1027 (1966). 

[7] B. M. Barker, M. S. Bathia and S. N. Gupta, Phys. Rev., 
158, 1498 (1967). 

[8] B. S. Dewitt, Phys. Rev. 160, 1113 (1967); ibid. 162, 
1195 (1967); ibid. 162, 1239 (1967). 

[9] K. A. Kazakov, Class. Quant. Grav. 18, 1039 (2001). 
[10] C. Itzykson and J. B. Zuber, Quantum Field Theory, 

(McGraw-Hill Book Co., Singapore, 1985) 
[11] V. N. Kondratyev, Phys. Rev. Lett. 88, 221101 (2002). 



■5 



[12] E. Kolb and M. Turner, The Early Universe (Addison- 
Wesley, 1988) 

[13] A "sufficiently high" photon density here means that the 
number of photons present in a cube of size A^ h , where 
A ch is the characteristic wavelength of the spectrum, is 
much larger than unity. Note that for a photon gas sub- 
ject to cosmological expansion, this number is a con- 
stant. Furthermore, for the particular case of the cosmic 



background radiation, this constant is much smaller than 
unity. 

[14] A. Hasegawa, Plasma Instabilities and Nonlinear Effects, 
Chap. 4.3. ( Springer- Verlag, Berlin, 1975). 

[15] A. Kallberg, G. Brodin, and M. Marklund, Class. Quan- 
tum Grav. 23, L7 (2006). 



